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Abstract. The influence of a thermodynamic constraint on the critical finite-size 
scaling behavior of three-dimensional Ising and XY models is analyzed by Monte- 
Carlo simulations. Within the Ising universality class constraints lead to Fisher 
renormalized critical exponents, which modify the asymptotic form of the scaling 
arguments of the universal finite-size scaling functions. Within the XY universality 
class constraints lead to very slowly decaying corrections inside the scaling argu- 
ments, which are governed by the specific heat exponent a. If the modification 
of the scaling arguments is properly taken into account in the scaling analysis of 
the data, finite-size scaling functions are obtained, which are independent of the 
constraint as anticipated by analytic theory. 

1 Introduction 

The theoretical investigation of classical spin systems has played a key role 
in the understanding of phase transitions, critical behavior, scaling, and uni- 
versality [|l]j^]. In particular, the classical Ising, the XY, and the Heisenberg 
model are the most relevant spin models in three dimensions. Each of these 
simple models represents a universality class which, apart from the spatial 
dimensionality and the range of the interactions, is characterized by the num- 
ber of components of the order parameter , e.g, the magnetization in the case 
of ferromagnetic models. Real systems, however, suffer from various kinds of 
imperfections, e.g., lattice defects, impurities, or vacancies. In an experiment, 
which is designed to probe critical behavior as a function of temperature, the 
presence of, say, impurities on the lattice constitutes a thermodynamic con- 
straint, because in a given sample the impurity concentration will remain 
constant during the temperature scans. According to the concepts of ther- 
modynamics the impurity concentration m can be written as the derivative of 
the grand canonical potential with respect to the chemical potential fa of the 
impurities, where other parameters like the temperature and the volume of 
the system are kept fixed. Now the question arises how the critical singular- 
ities in the grand canonical potential are affected when the thermodynamic 
ensemble is changed from 'fixed fa' to 'fixed fa\ where the location of the crit- 
ical temperature T c depends on the particular values of fa or fa, respectively. 
The answer to this question has been given a long time ago by Michael Fisher 
H . Provided, that the critical singularites have their usual form in the 'fixed 
[ii ensemble, then the constraint fa = const, amounts to a reparameteriza- 
tion of the reduced temperature t = {T — T c (ni))/T c (fa) of the constrained 
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system in terms of the reduced temperature r = (T — T c (/ij))/T c (/Xj) of the 
unconstrained system according to || 

t = ar + &r|Tr Q + ..., (1) 

where a and 6 are nonuniversal constants and the dots indicate higher order 
contributions. Apart from a linear term (|l|) contains a singular contribution 
which is characterized by the critical exponent 1 — a of the entropy density 
. Which of the two terms in (0) is the leading one for t,r — > depends 
on the sign of a. Within the Ising universality class in d = 3 dimensions 
a ~ 0.109 0] so that |r| ~ |£| !/ (!— Q ) to leading order and therefore the 
critical exponents f3 (order parameter), 7 (susceptibility), and v (correlation 
length) of the unconstrained system undergo 'Fisher renormalization' in the 
constrained system according to |3| 

P->F = p/0.-a), 7-Y = 7/(l-«). ^f' = *V(l-a). (2) 

The specific heat exponent a requires a more careful analysis, because the 
specific heat is the temperature derivative of the entropy which in addition 
to the 'renormalization' displayed in @ causes a sign change 

a — > a' = —a/(l — a). (3) 

Note that analytic background contributions to the entropy of the uncon- 
strained system become singular in the constrained system due to the singu- 
larity in the reparameterization given by (jl]) . 

Within the XY universality class in d — 3 the exponent a is negative , 
where probably the best current estimate a ~ —0.013 is obtained from an 
experiment on 4 He near the superfluid transition jj| . For negative a the linear 
term on the r.h.s. of (|l|) is the dominating one for t — * 0. However, the XY 
universality class a is so small, that in practice the singular term in (|l|) can 
never be neglected. Instead, the singular contribution to ([!]) gives rise to very 
slowly decaying correction terms which must not be confused with Wegner 
corrections to scaling . These correction terms have to be considered in any 
scaling analysis in order to obtain correct values for the critical exponents. 

If the system is finite, which is neccessarily the case for any Monte - Carlo 
simulation, all critical singularities are rounded, i.e., all quantities are analytic 
functions of the thermodynamic parameters |6| so that a thermal singularity 
as shown in (Q) does not occur. Critical finite-size rounding effects in, e.g., 
a cubic box L d are captured by universal finite-size scaling functions j(i|,|7| 
which restore all critical singularities in the limit L —> 00. Following the line 
of argument in || , (|l]) then has to be replaced by 

t = ar + T\rr a f(rL 1 ^), (4) 

where f(x) is the finite-size scaling function of the entropy density and x = 
tL x I v is a convenient choice of its scaling argument. For r — > at finite 
L the singular prefactor of f(x) in (Q) must be cancelled so that one has 
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f(x) = A\x\ a + . . . in the limit x — > 0, where A is a nonuniversal constant 
such that f(x)/A is a universal function of its argument. To leading order 
in r the reparameterization of the reduced temperature t of the constrained 
system is therefore linear in the reduced temperature r of the unconstrained 
system and one finds 

t = r(a + AL a/u ). (5) 

According to (^|) the finite-size scaling argument x in the constrained system 
is given by 

x = tL 1 ' v =tL 1 / u /(a + AL a / u ), (6) 

where the shape of the finite-size scaling functions is maintained [pi, i.e., 
the presence of the constraint only affects the form of the scaling argument 
x. For a > (^) asymptotically reduces to x — tL l l v /A for large L in 
accordance with Fisher renormalization (see (^)). For a < (j^) captures 
the aforementioned slowly decaying corrections to the asymptotic critical 
behavior in the XY universality class when a thermodynamic constraint is 
present. Note that A > for the Ising universality class and that A < for 
the XY universality class. 

In the remainder of this paper a simple spin model is introduced which 
can be efficiently simulated with existing Monte - Carlo algorithms both with 
and without constraints in three dimensions. For the Ising and the XY version 
of the model finite-size scaling according to (||) is tested for the modulus of 
the order parameter, the susceptibility, and the specific heat. 



2 Model and simulation method 

The model system which is investigated here can be described as an 0(N — 1) 
symmetric classical 'planar' ferromagnet in a transverse magnetic field. The 
model Hamiltonian reads 

u = -jj2J2 s ? s J- h J2 s ?> ( ? ) 

(ij) x=l i 

where (ij) denotes a nearest neighbor pair of spins on a simple cubic lattice 
in d = 3 dimensions. The lattice contains L lattice sites in each direction and 
in order to avoid surface effects periodic boundary conditions are applied. 
Each spin is a classical spin with N components Si — (S], Sf,..., S™) 
with the normalization \ Si\ — 1 for each lattice site i. The magnetic field h in 
only acts on the TV-components of the spins which are not coupled by the 
exchange interaction J. From the symmetry of the Hamiltonian it is obvious, 
that the model belongs to the 0(N — 1) universality class in d — 3, where 
nonuniversal quantities like the critical temperature T c — T c (h) depend on the 
strength of the transverse field h. Note that T c (h) is symmetric around h = 
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and decreases with increasing h, because the spins become more and more 
aligned with the A^-direction as ±/i is increased and due to the normalization 
condition the typical interaction energy between pairs of spins is decreased. 

The Hamiltonian given by (Q) defines the unconstrained model. The con- 
straint is imposed on the transverse magnetization M in the form 



The critical temperature of the constrained model is a symmetric and mono- 
tonically decreasing function of the prescribed transverse magnetization M. 
The transverse field h here plays the part of the chemical potential of impu- 
rities (see Sect. 1) and the transverse magnetization M accordingly plays the 
part of the impurity concentration rii. It is also possible to implement O(N) 
spin models with impurities or vacancies with diffusion in order to mimic 
the situation discussed in ||. However, the fact that the Hamiltonians given 
by (0) and (0) only require a single 'species' with a single coupling constant 
leads to some simplifications in the algorithms. Note that the symmetric con- 
strained model (M = 0) becomes equivalent to the symmetric unconstrained 
model (h — 0) for sufficiently large lattices. In particular, both versions of 
the symmetric model have the same T c . 

The Monte-Carlo algorithm is chosen as a hybrid scheme, where each 
hybrid Monte-Carlo step consists of 10 updates each of which can be one 
of the following: one Metropolis sweep of the whole lattice, one single clus- 
ter Wolff update ||, or one overrelaxation update of the whole lattice fLOf , 
where the latter can only be applied for N > 3. The Metropolis algorithm 
updates the lattice sequentially and works in the standard way for the un- 
constrained model. For the constrained model the constraint M = const, is 
observed locally by applying a Kawasaki update dynamics for the N th com- 
ponents of the spins. For each lattice site i a nearest neighbor site j is chosen 
randomly and a random amount of the TV*' 1 spin component is proposed for 
exchange such that + remains constant. Then new spin components 
(Sj, Sf, . . . , S^ 1 ) are proposed and the spin components (Sj,Sj, . . . , S^ -1 ) 
are adjusted according to the spin normalization condition = \Sj\ = 1. 
The local change AE of the configurational energy is calculated according 
to (||). According to detailed balance the proposed update is accepted with 
probability p(/3AE), where (3 — l/(ksT). For our simulation we have chosen 
p(x) = l/(exp(x) + 1). Note that all updates must be proposed such that the 
new spin at lattice site i is taken from the uniform distribution on the unit 
sphere in N dimensions. 




(8) 



where the Hamiltonian of the constrained model is given by 



N-l 




(9) 
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The Wolff algorithm also works the standard way || , except that only the 
first N — 1 components of the spins are used for the cluster growth, i.e., (0) 
and (|) are treated as planar ferromagnets. This means that a cluster update 
never changes the N th component of any spin so that the Wolff algorithm is 
nonergodic in this case. The cluster update is still a valid Monte-Carlo step in 
the sense that it fulfills detailed balance, however, in order to provide a valid 
Monte-Carlo algorithm it has to be used together with the Metropolis algo- 
rithm described above in a hybrid fashion. The use of Wolff updates allows 
us to take advantage of improved estimators @ for magnetic quantities. 

The overrelaxation part of the algorithm performs a microcanonical up- 
date of the configuration in the following way. The local configurational en- 
ergy has the functional form of a scalar product of the spins, where according 
to (0) and (||) only the first N — 1 components are involved. With respect to 
the sum of its nearest neighbor spins each spin has a transverse component 
in the (Sj,Sf, . . . , S^' 1 ) plane which does not enter the scalar product. The 
overrelaxation algorithm scans the lattice sequentially, determines this trans- 
verse component for each lattice site and flips its sign. This overrelaxation 
algorithm is similar to the one used in |l0| and it quite efficiently decorre- 
lates subsequent configurations over a wider range of temperatures around 
the critical point than the Wolff algorithm. However, overrelaxation can only 
be applied for N > 3. In the following only the cases N = 2 (transverse Ising) 
and N = 3 (transverse XY) are considered. 

In a typical hybrid Monte-Carlo step we use three Metropolis (M), seven 
single cluster Wolff ( C) updates for N — 2 and three Metropolis, five single 
cluster Wolff, and two overrelaxation updates (0) for N — 3 in the critical 
region of the models. The inividual updates are mixed automatically in the 
program so that the update sequences (M C C M C C M C C C) for N = 2 
and (M C C M O C M C C O) for N = 3 are generated as one hybrid 
Monte-Carlo step. The shift register generator R1279 given by the recursion 
relation X n = X n - p © X n - q for (p, q) = (1279, 1063) is used as the random 
number generator. Generators like this are known to cause systematic errors 
in combination with the Wolff algorithm However, for lags (p, q) as large 
as the ones used here these errors will be far smaller than typical statistical 
errors. They are further reduced by the hybrid nature of our algorithm due 
to the presence of several Metropolis updates in one hybrid Monte-Carlo step 
& 

The hybrid Monte-Carlo scheme described above is employed for lattice 
sizes L between L = 20 and L = 80. For each system size and temperature we 
perform at least 10 blocks of 10 3 hybrid steps for equilibration followed by 10 4 
hybrid steps for measurements. Each measurement block yields an estimate 
for all static quantities of interest and from these we obtain our final estimates 
and estimates of their statistical error following standard procedures. At the 
critical point (see below) two or three times as many updates have been 
performed. The integrated autocorrelation time of the hybrid algorithm is 
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determined by the autocorrelation function of the energy or, equivalently, 
the modulus of the order parameter, which yield the slowest modes for the 
Wolff algorithm. The autocorrelation times are generally rather short, at the 
critical point they range from about 5 hybrid Monte-Carlo steps for L = 
20 to about 10 hybrid Monte-Carlo steps for L = 80. The values for the 
equilibration and measurement periods given above thus translate to roughly 
100 and 1000 autocorrelation times, respectively. In order to obtain the best 
statistics for magnetic quantities a measurement is made after every hybrid 
Monte-Carlo step. All error bars quoted in the following correspond to one 
standard deviation. The simulations have been performed on the DEC alpha 
AXP workstation cluster at the Physics Department and on HP PJSC8000 
workstations at the Computer Center of the RWTH Aachen. 

3 Ising universality class 

For N = 2 (^) and (||) describe a classical Ising model in a (fixed) transverse 
field or with fixed transverse magnetization, respectively. In the following we 
will only consider the constrained model with the symmetric constraint M = 
and with the constraint m = M/L 3 — 1/ y2. The symmetrically constrained 
model does not show Fisher renormalization ||] and we therefore use this case 
for tests of the algorithm and for the production of data representative of the 
Ising universality class in d = 3. The constraint m = const. ^ breaks 
the £j — * —Sf symmetry of the model and Fisher renormatization should 
become visible within a certain temperature window around T c = T c (m). The 
width of this window is of course a nonuniversal property of the model and 
in particular one expects this window to widen as m is increased. Due to the 
spin normalization condition m cannot exceed unity and one therefore also 
expects, that critical behavior becomes very difficult to resolve numerically 
if m is too close to its maximum value. Therefore, m = l/v2 is chosen as a 
compromise between good resolution in the critical regime and a prominent 
Fisher renormalization effect. 

The critical temperatures T c (m = 0) and T c (m — l/v2) are determined 
from temperature scans of the Binder cumulant ratio according to standard 
procedures H|. We obtain the following reduced critical coupling constants 
K c (m) = J/k B T c (m): 

K c (0) = 0.41638 ±0.00005 and K c (l/V2) = 0.6371 ± 0.0001. (10) 

The corresponding estimates for the Binder cumulant ratio obtained for 
m = and m = 1/V2 agree with previous estimates obtained for the Ising 
universality class within two standard deviations jlj|, where for the latter 
choice of m Wegner corrections to scaling are considerable and must be sub- 
tracted in order to obtain a reliable estimate. In order to obtain an estimate 
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for the exponent v which enters the finite-size scaling argument according to 
(^|) the cumulant 

has been measured, where <fi = L~ 3 ^\ Sj is the order parameter. At the 
critical temperature T c (m) the scaling behavior X ~ x/t is expected (see 
(0)). Corresponding numerical results for m = and m = 1/V2 are displayed 
in Fig|l] on a double logarithmic scale. The data are compatible with simple 

1000 f 1 




100 



Fig. 1. Cumulant X at the critical point for m = (x) and m = l/v2 (+). The 
solid and dashed lines display power law fits to the data for 30 < L < 70 for m — 
and m = 1/ \/2, respectively 

power laws, where the exponents v — 0.622 ± 0.005 (m = 0) and v' = 
0.714±0.004 (m = 1/ V2) have been obtained. Compared to the best currently 
known estimate v ~ 0.630 Q the above estimate is too small and only agrees 
with the theoretical value within two standard deviations. A more thorough 
analysis shows that the discrepancy can be explained by a mismatch of the 
order 5 x 10~ 5 between the actual critical temperature and the estimate 
used here (see (fL0|)), which on the other hand is of the same magnitude as 
the statistical error of K c (0). The agreement between the above estimate for 
v 1 = v/(l — a) and the theoretical value v' ~ 0.708 H is better, however, it 
may again be affected by a mismatch between the actual value of T c (l/y/2) 
and the estimate used here. If the literature values for v and a are substituted 
in (^|), where a and A are used as fit parameters, a/ A ~ 0.1 is obtained which 
is small enough to be ignored in the scaling analysis (see below). 
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The finite-size scaling analysis has been performed for several thermody- 
namic quantities, in particular, the average modulus of the order parameter 
the susceptibilities 

X+ ^^(0 2 ), x- = ^ «0 2 } - (I0l> 2 ) , (12) 

and the specific heat C. Data will only be shown for (|</>|), X-j an d C, because 
the finite-size scaling functions for (\<fi\) and \+ are very similar. According 
to finite-size scaling theory it must be possible to callapse the data for all 
m onto one and the same curve, where two nonuniversal scaling factors are 
required for each quantity One scaling factor adjusts the magnitude of the 
scaling argument x (see (pi)), the other adusts the absolute normalization of 
the quantity. Note that the former saling factor must be the same for all 
quantities. For m = the scaling argument x — tL 1 ^ is used, whereas for 
m = 1/V2 the choice x = tL( 1 ~ a ^ l/ /A has been made, where A ~ 1.1 and 
the coefficient a in (6) has been neglected. The exponents v and a are taken 
from the literature [4]. The scaling plot of (\4>\) is shown in Fig.pl where the 




Fig. 2. Scaling plot of (\<j>\) for L = 30, 40, 50, and 60 for m = and m = l/y/2. 
The reduced temperature t has been varied between —0.007 and 0.003. Statistical 
errors are much smaller than the symbol sizes 



abolute normalization of the data for m — l/\/2 can be adjusted to the m = 
data by a scale factor of ~ 0.7 as one would expect from simple mean field 
arguments. As shown in Fig.^, the espected data collapse can be reproduced 
rather well, where the literature value for the exponent (3/v — 0.5168 j| 
has been used. The same holds for the susceptibility X- which is displayed 
in Fig||, where j/is = 2 — 77 = 1.967 is also taken form B. The absolute 
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magnitudes of x-fL 1 ^ for m — and m — l/\/2 are different by a factor 
of about 0.5 which is in accordance with simple mean field arguments. Note 
that the scaling function of x- has a maximum for x < ]16[ . 
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Fig. 3. Scaling plot of X - for L = 30, 40, 50, and 60 for m = and m = 1/V2. 
The reduced temperature t has been varied between —0.007 and 0.003. Statistical 
errors are much smaller than the symbol sizes 



The specific heat C requires a somewhat different treatment due to the 
fact that unlike the other quantities presented so far the specific heat requires 
an additive renormalization within renormalized field theory pq |. For the 
data analysis this means that scaling can only be obtained after a suitable 
subtraction is applied to the specific heat. One option to obtain scaling is 
to subtract the bulk specific heat C°(t) at a reference reduced temperature 
t = to which are given by 

C?(*) = —\t\- a + B and t = {L/e±T 1/v , (13) 
a 

where A± and B are nonuniversal constants and £± is the amplitude of the 
correlation length. The index ± refers to temperatures above or below T c (m), 
respectively. The reference reduced temperature chosen in ( |l3| ) is positive 
and therefore only A + and are needed. Specifically, the choice v = 0.630, 
a = 0.109 §, A+ = 0.1552, B = -1.697, and £° = 0.495 @ guarantee 
scaling of the relative specific heat AC = C — CP(to) for the Ising model in 
d = 3. Note that £9 is measured in units of the lattice constant. For the data 
to be analyzed here the subtraction defined by ( Jl3| ) is only valid for the case 
m = 0, where AC scales as L a l v . For m = Fisher renormalization 
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according to (|l|) must be applied to C®(t) in order to obtain the correct form 
C™ (t) of the subtraction. The result is 

CTW = ^ItH^-^ + B' and t = (i/rf^"" 1 ^ , (14) 
a 

where A' + = -0.1728, B' = 1.598, and £° = 0.495 is not changed. The 
scaling factor L a ^ , which usually governs the finite-size scaling of the specific 
heat, is cancelled here, i.e., one expects data collapse for AC°/L a / u and 
AC' m = C — C™(to) up to an overall scale factor of about 0.5. The result of 
the data analysis is shown in Fig.|[ The data collapse reasonably well onto 
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Fig. 4. Relative specific heats AC°/L a/ " for m = and AC m for m = 1/V2 for 
L — 30, 40, 50, and 60. The reduced temperature t has been varied between —0.007 
and 0.003 



a single curve except near the maximum of the scaling function, where also 
the scatter of the individual data is substantial due to a few bad samples for 
L = 50 and 60. However, there are also systematic deviations from scaling in 
the data, because the maximum in the scaling functions for the m — data 
is more pronounced than in the m = l/y/2 data. These deviations could be 
due to enhanced Wegner corrections to scaling for m = l/v2 as compared 
to m = 0. 



4 XY universality class 



For N — 3 (0) and (|9|) describe a classical XY model in a (fixed) transverse 
field or with fixed transverse magnetization, respectively. As for the case 
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N = 2 we will only consider the constrained model with the symmetric 
constraint m = and with the constraint m = 1/V2 in the following. The 
symmetrically constrained model is again used for algorithmic tests and data 
production for the XY universality class. The nonsymmetric constrained XY 
model does not show Fisher renormalization, however, according to ([|) very 
slowly decaying corrections to the asymptotic critical behavior are expected, 
which will be discussed in the following. First, the cumulant X defined by 
( p"T| ) is evaluated at the critical point, which is given by the reduced coupling 
constants 



K c (0) = 0.6444 ± 0.0001 and K c (l/V2) = 1.1126 ± 0.0003, (15) 
respectively. The result for X is displayed in Fig.||. For m = the data can 
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Fig. 5. Cumulant X at the critical point for m = (x) and m = l/y/2 (+). The 
solid and dashed lines display fits to the data for 30 < L < 80 for m = and 
m = 1/V2, respectively (see main text) 



be fitted by a power law ~ L x l v , where v = 0.678 ± 0.008 is obtained which 
agrees with the best current estimate v = 0.671 Q. For m — the data 

can also be fitted by a power law, however, the resulting exponent v only 
has the meaning of an effective exponent which does not fit into the XY 
universality class. As shown in Fig.|5| the expression x/t according to (^|) also 
yields a very good representation of the data where the parameter b in FigJ^ 
is given by b = A/a = —0.941. The exponents used in the fit (XY universality 
class) are taken from Q . The value of the Binder cumulant found here agrees 
with results reported in the literature for the standard (plane rotator) XY 
model ]l7| ], however, Wegner corrections to scaling become quite substantial 
for m = 1/V2. 
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In the following scaling analysis the finite-size scaling argument for the 
case to = takes its standard form x = tL 1 ^ and for to = l/y/2 the 
combination x — tL 1 ^ /(l + bL a / v ) for b = —0.941 takes care of the slowly 
decaying correction terms to the asymptotic critical behavior caused by the 
very small and negative value of a in the XY universality class. As in Sect. 3 
we consider (|(/>|), X+i X—i an d the specific heat C in the scaling analysis. 
The scaling functions for (\<p\) and x+ again look very similar so that we 
do not reproduce scaling plots for x+ here. The result for (|^|) is shown in 
Fig|| where the order parameter 4> = L~ 3 Y^ii^l , Sf) has two components 
here. The data collapse very well onto a single curve. Note that the absolute 




Fig. 6. Scaling plot of (\(j>\) for L = 30, 40, 50, and 60 for m = and m = l/y/2. 
The reduced temperature t has been varied between —0.005 and 0.005. Statistical 
errors are much smaller than the symbol sizes 



magnitudes of (\(f>\) for m = and to = l/y/2 are again related by a factor 
of ~ 0.7 as suggested by mean-field arguments. The values for the critical 
exponents v = 0.671 and /3 = 0.347 are taken from the literature j|. The 
susceptibility X- can be treated essentially as described in Sect. 3, where the 
exponent 7/1/ = 2 — r\ = 1.965 is taken from [0. The result of the scaling 
analysis is displayed in Fig.f7[ The data do not collapse as well as in Fig.[| 
Especially near the maximum of the scaling function the scatter of the data is 
substantially larger than in FigM Slight systematic deviations from scaling 
for to = 1/V2 are observed which may again be due to enhanced Wegner 
corrections to scaling as compared to to = 0. Note that contrary to Fig. [| the 
scaling function of x~ has a maximum for x > in the XY universality class. 

The specific heat C of the XY model also requires a subtraction before 
scaling is obtained [p|. The subtraction C°(io) is again used in the form 
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Fig. 7. Scaling plot of X - for L = 30, 40, 50, and 60 for m = and m = 1/V~2. 
The reduced temperature t has been varied between —0.005 and 0.005. Statistical 
errors are much smaller than the symbol sizes 

given by (§J|), where A + = 0.42, a = -0.013, B = -A+/a, and £° = 1.0 are 
used here which differ somewhat from the choices made for the standard XY 
model in fTs)] . It turns out, that the quality of the data collapse for the relative 
specific heat AC = C — (to) for m = is rather insensitive to the choice 
of Q. The form of the subtraction C b m (i ) for m = l/y/2 requires a little 
analysis in order to include the slowly varying corrections to the asymptotic 
behavior coming from (|l|). One obtains the approximate form 

<o = L- 1/l/ (l + 6£ Q/ "),(16) 

where b = —0.941 (see FigJ|), A+ = 0.42 as before, and c ~ 2.0 for optimal 
data collapse. The resulting scaling plot of AC and AC' m = C - C b m (i ) 
is displayed in FigJ|. The overall shape of the scaling function is similar to 
the one shown in Fig.^J. However, the scatter of the data near the maximum 
is so strong, that data collapse cannot be obtained in this region. In part 
this deficiency in the data may be due to the presence of 'bad' samples, e.g, 
for L = 50 and m = at t = -0.003 and t = -0.002 and for L = 60 
and m = l/y/2 at t = —0.005 and t = —0.003. Apart from that deviations 
from scaling as in Figji] may be present which are due to an enhancement of 
Wegner corrections to scaling for m — l/y/2 as compared to the case m = 0. 
However, Figs.^ - ^ confirm, that the choice of the scaling variable x given by 
(||) captures the slowly decaying correction term inside the scaling argument 
in an appropiate way and that furthermore the finite-size scaling behavior of 
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Fig. 8. Relative specific heats AC (0 ' m) / [L a/u /(l + bL a/u )\ with b = for m = 
and 6 = —0.941 for m = 1/V2 for L = 30, 40, 50, and 60. The reduced temperature 
t has been varied between —0.005 and 0.005 



constrained models in the Ising and the XY universality class can be treated 
on the same footing. 

5 Summary and conclusions 

The influence of constraints on the critical finite-size scaling behavior of Ising 
and XY models has been investigated by Monte-Carlo simulations of O(N-l) 
planar ferromagncts with fixed transverse magnetization. The theoretical idea 
that only the form of the scaling argument is modified, whereas the shape 
of the universal scaling functions remains unchanged is verified within the 
statistical uncertainty of the data for the modulus of the order parameter, 
the susceptibilites x+ (not shown) and and the specific heat. The form of 
the scaling argument used here allows to deal with critical finite-size effects 
in constrained Ising and XY models on the same footing, where constrained 
Heisenberg models can be included as well. Within the Ising universality class 
the finite-size behavior is consistent with the Fisher renormalization of crit- 
ical exponents. In the XY universality class slowly decaying corrections to 
the asymptotic critical behavior are generated which are captured systemat- 
ically by the analytic form of the scaling argument. The treatment of these 
corrections within the XY universality class may serve as a paradigm for 
the finite-size scaling analysis of dynamic quantities, where the constraints 
imposed here reappear as conserved quantities which are statically or dy- 
namically coupled to the order parameter. These corrections may also be 
important for the interpretation of spin dynamics data for planar ferromag- 
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nets, where the energy of the system is conserved during the simulated time 
evolution of spin models similar to the ones investigated here. 
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